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Abstract 

We propose to test the CP symmetry in the reactions Z — > fi + fi~j 
and Z — > T + r^'~f. The experimental analysis of angular correlations 
allows to determine a set of effective couplings: the electric and 
weak dipole moments of the muon and the tau lepton, and in par- 
ticular chirality conserving 4-particle couplings, all of which can be 
induced by CP-violation in renormalizable theories of electroweak 
interactions beyond the Standard Model. We update an indirect 
bound on the weak dipole moment of the muon. 
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1 Introduction 



The search for CP violation in Z boson decays provides an interesting tool to 
investigate possible new physics beyond the Standard Model. Several proposals 
have been made how to look for CP violation in processes such as Z — > £ + £~, 
£ + £~7, qq, qq7, qqG, where t = e, /x, r leptons, q = u,d,s,c,b quarks, and G 
= gluon (cf.[l]-[12] and references therein). CP-violating effects in the above 
reactions can be parametrized by form factors, among them the well known 
electric and weak dipole moments of quarks and leptons, or one can use the 
effective Lagrangian approach. Experimental investigations of CP-odd effects 
in Z decays have been performed resulting in upper limits on the weak dipole 
moment of the r lepton [13]— [15] and on a CP violation parameter in Z — > bbG 
[16]. An "indirect" limit on the r electric dipole moment was given in [17]. 

Today, a large number of Z bosons has been produced in electron-positron 
collisions both at LEP (CERN) and at SLC (SLAC). Thus investigations of 
relatively rare events like radiative three-body decays are possible. In this 
article we investigate in detail the unpolarized electron-positron annihilation 
to a muon or tau lepton pair and a hard photon with respect to CP violation: 

e+(p + ) + e"(p_) - Z - t{k l )+ e-(k t )+ 7 (g), (£ = »,t). (1) 

Within the framework of the Standard Model (SM), CP violation in these 
leptonic reactions occurs only at higher orders in perturbation theory, leading 
to unobservably small effects. However, interactions beyond the SM such as 
the exchange of virtual photinos or other neutralinos, neutral Higgs bosons, 
excited leptons, or leptoquarks, can induce CP-violating contributions to the 
interaction of leptons and gauge bosons. A calculation of the dipole form 
factors of the r lepton in several extensions of the SM can be found in [18]. 
Most of these models have in common that they lead to CP-violating effective 
couplings, typically generated at one-loop level, which are are not flavour- 
universal but can grow substantially with the mass of the quarks and leptons 
involved. 

In the following section a model-independent parametrization of CP-odd 
couplings relevant for the reactions (1) is discussed. Then we propose some 
easily measurable CP-odd momentum correlations which are optimized to ex- 
tract information on these couplings from the data. In sections 4 and 5 we 
discuss in some detail the level of accuracy of CP tests obtainable in Z decays 
to /i + /i~7 and r + r~7, respectively. 
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2 Effective Lagrangian approach 



In order to describe in a model-independent way possible CP-violating effects 
we use an effective Lagrangian with CP-odd operators of mass dimensions 
d > 4. Restricting ourselves to d < 6 (after symmetry breaking!), the operators 
which are relevant for the reactions considered here are given by (c.f. [2, 6]) 



£cp(X) 



= E{ 

e 

l -dJ £{x)a^ lb £{x) [d^A v (x) - d„A^x)\ 

l -d* i(x)a^ 5 e(x) \d,z u (x) - d u z,(x)] 



(2) 



+ 



f ve e(x)Y£(x) + f M l{x)Yld{x) Z»{x) [d,A u (x) - d v A^x)]}, 



where dj, df are the electric dipole moment (EDM) and weak dipole moment 
(WDM) coupling parameters of the leptons £ — /i,t and fve, f At are real cou- 
pling constants corresponding to CP-violating Z-photon-lepton 4-particle ver- 
tices. The photon and Z-boson fields are denoted by A^ and Z M , respectively. 
The authors of Ref. [11] have given a form factor decomposition for reactions 
of the type (1). In such a framework our couplings in (2) correspond in essence 
to the CP-violating form factors with minimal momentum dependence. 

Following previous analyses [2, 8] we introduce the dimensionless parameters 
fve, fAt, dj, df defined by 



fve,M — 
d]> z = 



e 2 ■ Q e c 

- — 7. 7. 2 iV£,A£ 

sin u w cos v\v m z 



2.62 • 10" 
GeV 2 



-fve, A£i 



e • Qe 



sin 6 W cos 0w m z 



d 



<7,Z 



= 5.15- 10~ 16 ecm-^, 



(3) 
(4) 



where £ — /i,t and Qe — — 1, and e 2 / (47r) is the fine structure constant. We use 
m z = 91.19 GeV and sin 2 6* w = 0.23 for the numerics. We work to leading order 
in the CP-violating couplings of C C p and in the couplings of the SM. Then 
the parameters in the effective Lagrangian are identical to the corresponding 
form factors, and imaginary parts of the form factors are absent. This point 
was discussed in detail in [12], and it was shown that in this order one has 



Re d 



J7,Z 



Im dj'^s) 



J7,z 
a e ' 

0. 



(5) 
(6) 
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Here dj ,z (s) are the EDM and WDM form factors measurable in the reaction 
e + e~ — > at cm. energy y/s. We note that in a pure form factor approach 
the form factors relevant for e + e~ — ■> and e + e~ — > £ + £~ , ~f are not directly 
related. Such a relation is only achieved after making an assumption for the off- 
shell dependence of the dipole moment form factors. On the other hand, in the 
effective Lagrangian approach there is no problem to relate the parameters of 
the above two reactions in leading order. Also in higher orders it is in principle 
straightforward to work out this relation (cf. [12] for further remarks). 

The next step is to calculate the differential cross sections for e + e~ — > 
and e + e~ — > i + i~ r y taking into account the SM diagrams (one is shown in 
Fig. la) and the diagrams induced by the couplings of £cp (Figs- lb,c,d). 
Neglecting radiative decays with more than one 7 the sum of the two channels 
above gives the inclusive width 

r(z - e+rx) = r SM (z - e + rx) + Ar CP , (7) 

where r SM is the SM width and Ar CP is of second order in the CP-violating 
couplings of Cqp- We use this calculation to make an indirect estimate of the 
order of magnitude which is allowed for the CP-violating parameters in Cqp 
(cf. (2)). The measured widths and T TT which include all radiative events 1 
agree well with the predictions from the SM (see Tab. 1). 





r 

1 exp 




jsriv + 5Tl u 


n + n X 

T+T-X 


83.79 ± 0.22 MeV 
83.72 ± 0.26 MeV 


83.97 ± 0.01 ± 0.06 MeV 
83.97 ± 0.01 ± 0.06 MeV 


0.23 MeV 
0.27 MeV 



Table 1 Experimental and theoretical (SM) values of the total decay widths 
for Z -> fi+fi-X and Z -> t + t~X (cf. [19, 20]). 



We define values Ar cxp for the "allowed" extra contribution to the width at A 
standard deviations as 

at££- = r cxP - r SM + a • y^rLp + sr^. (8) 

In principle a combined estimate of the four CP-violation parameters of (2) 
for each lepton £ is possible by ascribing to them any deviation from the width 
as predicted by the SM, i.e. by setting Ar CP < Ar cxp . However, this would 

1 We thank N. Wermes and M. Wunsch for providing information concerning this point. 
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ignore the possibility of other non-standard but CP-conserving effects to the 
widths which can lead to an increase or decrease of the widths as compared 
to the SM. Thus such a procedure only gives an order of magnitude estimate 
of the values allowed for the couplings in Ccp- We give a simplified estimate 
by ignoring interference terms induced by the new couplings, i.e. we assume 
that only fy e and f^i or dj or of contribute in the width T(Z — > £ + £~X). As 
compared to the other couplings, the contribution of the electric dipole moment 
dj is suppressed by a factor of a ew in the electroweak perturbation series. Thus 
the leading terms in APcp, calculated without phase space cuts and neglecting 
the lepton masses, are easily obtained from Tab. 1 of Ref. [8] as follows 

Ar CP (df) = ( — M! ) -0.24MeV, (9) 

\ 10 Lt ■ ecm J 

Ar CP (WA£) = + 0.042 MeV. (10) 

Now we can derive A s.d. limits on df and fa, f^e by requiring the theoretical 
widths (9,10) to be smaller or equal AT%£- (8): 

Ar CP (rff) < Art p d -, (11) 

Ar CP (/ w ,/ A ,) < ArX; (12) 

We see from Tab. 1 that the differences of the central values, r exp — r SM are 
negative. Correspondingly the 1 s.d. values AT^" are very small. Limits on 
df, fve an( i fpd obtained from (11, 12) at 1 s.d. are then not to be trusted. As 
explained above, we should consider (11, 12) not as strict inequalities but only 
as indications for the order of magnitude which is possible for the CP-violating 
couplings. We thus prefer to quote the 2 s.d. limits below, since then Ar^f' is 
at least of the order of the combined experimental and theoretical errors given 
in the last column of Tab. 1. 

In this way we get from (11) for the WDM of the muon the 2 s.d. upper 
limit 

\dl\ < 0.021 or \dl\ < 1.1 • 10" 17 ecm, (13) 
while fyg and f M are restricted from (12) at the 2 s.d. level by 

/v„ + /X. < ( 2 - 6 ) 2 > ( 14 ) 

fir + fl < (2-6) 2 . 
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This type of reasoning to obtain limits for the dipole moments was frequently 
used (cf. e.g. [4, 8]). The estimate for dP in (13) updates the one given in [4]. 

For the WDM of the r lepton it is not necessary to rely on indirect limits 
since experimental limits exist, obtained from direct CP violation searches. We 
quote the latest result for |Red^| which is related to our rff as explained above: 
The combination of all LEP results [21] gives the 95% c.l. upper limits 

\df\< 0.007 or |Rec^| < 3.6 • 10~ 18 ecm. (15) 

Estimates of the EDMs which contribute at tree level to radiative events 
only (see Fig. 1) can be given in a similar way by exploiting the experimentally 
allowed range between the SM prediction and the measured values of T(Z — > 
£ + £~ , ~f). An upper limit was given for the r EDM in [22]. Using the Z decay 
data collected up to now it will be possible to update that value. The EDM of 
the muon is known to be smaller than 10~ 18 ecm, or \d?\ < 0.0019 [20]. (Again, 
we can identify the on-shell EDM with the corresponding coupling parameter 
in (2) for our purposes.) The momentum correlations proposed in the following 
section are blind to the EDM, therefore our analysis of Z — > fi + n~ / ~f aims at 
only three still unknown parameters /v M , /a^ and d 7 . 

The calculation of the differential cross section shows that it is convenient 
to choose instead of the parameters fye and f^e the linear combinations 

fa = 9ve ■ Ia£ - 9ki • fve, (16) 
fa = 9ve ■ fve — 9 hi • fhe, 

where gye = — 2Q e sin 2 9w and g^e = T 3 £ = —1/2 are the usual neutral 
current coupling constants of the leptons £ = /i,t. In the zero-mass limit fu 
induces an angular distribution of the £ + £~ , ~f momenta which is measurable with 
tensor observables, while fa is related to vector observables. Thus, with ap- 
propriate observables these couplings can easily be determined independently. 

3 Observables 

In the following we consider the case of unpolarized e + e~ beams in the reactions 
(1). Then the initial state is described in the cm. system by a CP-invariant 
density matrix (CP tests with longitudinally polarized beams have been dis- 
cussed in [18, 23]). Furthermore we assume that the experimental phase space 
cuts and the momentum reconstruction methods do not introduce a CP bias. 
Then the mean values of CP-odd momentum correlations trace the couplings 
defined in £cp- 
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The task is now to construct a set of momentum correlations whose distribu- 
tions are sensitive to the couplings defined in £ CP with a reasonable statistical 
significance. For simplicity, we first consider each of the parameters fu, fu, 
df and dj separately by setting the other three to zero. Then the linear con- 
tribution to the mean value of a given observable, together with the width of 
its distribution, yields an estimate of the statistical accuracy with which this 
coupling can be measured. 

The kinematics of the reaction Z — > /i + /i~7 is determined by the fj, + , fx~ mo- 
menta k ± in the Z rest frame which is to a good approximation the laboratory 
system at LEP and SLC. In the case of Z — > t + t~7 with subsequent one- 
prong r decays which are the most frequent ones (about 85 %), the measurable 
quantities are primarily the momenta of the charged decay products, again 
denoted by k±, and the photon momentum q. Using these momenta one can 
construct many CP-odd observables. We have calculated the signal-to-noise 
ratios of some basic correlations and found the following observables especially 
useful: 

T — (k+ — k_) • P+ (k+ x k_) ■ P+ , (17) 
V = (k + xL).p + . 

Here hats denote unit momenta, and p + is the initial positron momentum vec- 
tor. With respect to the final state momenta T is a component of a tensor 
and V a component of a vector observable. These correlations involve neither 
particle energies nor the photon momentum. For Z — > r + r _ 7 we addition- 
ally studied the observables T* and V* which are defined as above but with 
k ± taken in the rest system of the t + t~ pair. The r momenta usually can- 
not be fully reconstructed from the observed decay products and the r decay 
vertex information. But the observation of the photon should be sufficient to 
reconstruct in each event the r + r~ center of mass frame. 

Let us denote the coupling parameters fie, he, df , dj (cf. (2-4, 16)) by 
gie, ...,gu- The mean values of CP-odd observables such as T and V defined 
above are in general not purely linear functions of the gn since the normalization 
introduces a denominator with quadratic contributions in the gus. The second 
order terms certainly are negligible if the couplings are small. For practical 
purposes it might be useful to consider observables which get linear contribu- 
tions only. These can easily be obtained by multiplying any of the momentum 
correlations discussed in this article with the measured width T(Z — > 
for the phase space region, i.e. the cuts considered. In this way one con- 
structs observables whose mean values are strictly linear in the couplings to be 
measured: 
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r(Z-*+r 7 ) • (O) = j^ai-gu. (18) 

1=1 

The Oj are proportionality constants which can be calculated numerically, and 
O denotes any CP-odd observable which does not depend on the unknown 
couplings. 

In order to get an impression about the quality of the observables intro- 
duced above we also investigate for the couplings ga the corresponding optimal 
observables which have the greatest possible statistical signal-to-noise ratio. 
For a single coupling this method was introduced in [24, 25]. The generaliza- 
tion to an arbitrary number of new couplings was given in [26] (cf. also [27]). 
The optimal observables are constructed as follows: Let us write the differential 
cross section as 

4 4 

da(<f>) = da SM (<f>) + ]T da^i^gu + ]T dag } (0)^. (19) 

i=l i,j=l 

Here stands collectively for the phase space variables, da S M denotes the CP- 
conserving SM part of the differential cross section, and daf^ are the contribu- 
tions due to the interference of the CP-odd interactions of £cp with the SM 
amplitudes. Then the ratios 

°<<*> = Sf (1 - ' - 4) (20) 

are the observables with optimal statistical sensitivity for a joint measurement 
of the couplings assuming that the couplings which shall be measured are small. 
The observables Oi (20) are optimal for arbitrary phase space cuts. Their 
sensitivity, however, depends on the cuts. 

For the reaction Z — * /i + /i - 7 one can easily measure the momenta of all 
final state particles. Let q be the 7 momentum and kp = k + ,k^ = k_ the 
momenta of /i±. Thus in this case the phase space variable in (20) is given 
by 

0=(k + ,k_,q). (21) 

Due to the energy-momentum conservation constraint the phase space is 5- 
dimensional. Identification of a photon requires some isolation cuts. For our 
analysis we assume covariant cuts of the form 

M±M! > y, (22) 



rn 



z 
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where fc(i), k@) are any two different four momenta from the set k + ,k-, q, and 
y is chosen to be 0.03. 

For the reaction Z — > t + t~ 7 with subsequent 1-prong r decays 

e+(p+) + e"(p_) - Z - r+(fc f ) + r"(A; T ) + 7 (g), (23) 



r + (/c f ) -> a + (£;+)+ neutrals, (24) 
T~(k T ) — > b^(k^) + neutrals, (a,b — tt, p,e, fx) 

we assume integration over the phase space variables of the neutrals. We are 
then left ideally with the information from the momenta of r + , r~, 7 , a + , 6~ 
and thus here we can identify 

= (k ? ,k r ,k + ,k_,q). (25) 

Of course, these momenta are not totally independent. Energy-momentum 
conservation in production and in the decays imposes well known constraints. 
Depending on the decay channels (24) the phase space is 9 to 11 dimensional. 
For our analysis below we assume again the isolation cuts (22) for the momenta 
of the 7 and the charged decay products a + ,b~. In practice the complete in- 
formation on the r ± momenta may not always be available. We will, therefore, 
discuss the optimal observables for the ideal case where the r ± momenta can 
be reconstructed. This then gives the ideal sensitivities. On the other hand, 
for the simple observables T, V of (17) and their analogues T* , V* constructed 
in the t + t~ cm. system one does not need knowledge of the individual r + and 
r~ momenta. Obviously, with partial information on the r ± momenta one can 
then reach sensitivities in between. 

From (20) we get as many optimal observables as we have parameters. Of 
course, in general the expectation values (Oi) will get linear contributions from 
all four parameters in leading order. We have investigated the matrix describing 
this dependence of the mean values of observables on the parameters. For the 
optimal observables this is in leading order identical to the correlation matrix 
of the Oi (cf.[26]). We found that these correlations are either negligible, or 
the couplings contributing are bounded by independent experiments. In the 
following we therefore confine ourselves to discuss the diagonal elements of the 
covariance matrix. We calculated numerically the mean values and the widths 
of the Oi and got in this way for a given number of events the statistical 
accuracy with which each of the unknown parameters can be measured ideally, 
assuming that all other new couplings are zero. For a numerical evaluation it is 
not necessary to have the Oi(<f>) in explicit form. We use a Monte Carlo event 
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generator 2 for the reactions (1) which takes into account the effects of £cp- 
Since the optimal observables are composed from terms of the differential cross 
section, one can easily generate them numerically with a Monte Carlo method. 

The knowledge of the optimal observables is interesting by itself, and useful 
for deriving simpler correlations which are approximately optimal. In most 
cases we found such observables which presumably are preferable for an exper- 
imental analysis. 



4 Z decays to n + n "7 

The LEP experiments have collected several thousand events of the reaction 
e + e~ — > /i + /i"7. For the CP tests proposed here it is sufficient to measure the 
muon directions of flight k±. Using these one can calculate the mean values 
of the observables T and V defined in (17). The tensor correlation T gets a 
contribution only from allowing a clean determination of this parameter. 
On the other hand, the mean value of the vector observable V depends both on 
f 2f M and a^. As expected, the vector correlation is slightly less sensitive than the 
tensor observable, due to the smallness of the vector polarization of the Z boson. 
The contribution of the weak dipole moment to V is negligible: in order to cause 
a measurable effect it would have to be three orders of magnitude bigger than 
the indirect upper limit (13). We conclude that f 2fM can be measured practically 
without contamination by other CP-violating sources using the observable V 
of (17). 

We also investigated the correlation of a measurement of our parameters 
with optimal observables. Since the EDM coupling does not contribute to 
CP-odd terms in leading order (cf.[2]), we have here three parameters f ifl , f 2fl , 

(denoted by g ltl , (? 3jU ) and three observables Oj (cf. (20)). Working always 
in leading order in the g ifl , we have 

(Oi> = Ecy ffj> (26) 

3=1 

with 

Cij = (OiOj)sM- 

2 Monte Carlo event generators for the reactions e + e 
with subsequent r decays are available from the authors. 



[i + H 7 and e + e 



(27) 



r+r <y 
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We calculated this covariance matrix. In Tab. 2 we give the cu (i— 1,2,3) and 
the rescaled matrix 

(O l O,) SM / v ^-c-, (28) 

where the diagonal elements are equal to 1. From these numbers we can easily 
derive parameters g[ {% = 1, 2, 3) where the corresponding optimal observables 
0\ have a diagonal covariance matrix c^- = (O^O'^sm- These parameters can 
be chosen as 



g[ — fin-, 

92 = /2„ - 1.1 • 10- 4 <5, (29) 
9' 3 = dl 

Recall that for optimal observables the covariance matrix V(g') of the estimated 
couplings is the inverse of the number of events times the covariance matrix of 
the observables (cf. [26, 27]) 

V(g>) = (No')- 1 - (30) 

Thus a diagonal d ensures a diagonal V(g'). We conclude from (29) that also 
for optimal observables the correlations are negligible. For a general discus- 
sion of diagonalization procedures in the estimation of more than one coupling 
parameter we refer to [27]. 

From the above we find that we have to consider for any of our observables 
O = T, V or an optimal one only a dependence on one coupling parameter. 
Calling it generically g we get in linear approximation 



(0) = c-g. (31) 

The 1 s.d. accuracy obtainable for a determination of g by a measurement 
of O can be estimated as 

S g =^, (32) 



where (C 2 )sm is the variance of O, calculated in the SM, and N is the number 
of events within the cuts considered. 

In Tab. 3 we summarize the statistical accuracies Sg of (32) with which 
the WDM of the muon, f ltl and / 2(U can be measured assuming a sample of 
iV = 10 4 events within the cut (22). The numbers in Tab. 3 represent 
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the minimal sizes of the couplings which lead to visible effects at f s.d., or, 
if no deviation from the CP symmetry is found, Sg is the 1 s.d. upper limit 
on the coupling parameter. The upper limits on / 2/1 correspond to two 
stripes in the /v M , /a^ parameter space. These are depicted in Fig. 2, together 
with the indirect bound (14) assuming zero mean values and a phase space cut 
parameter y = 0.03. The sensitivities of the optimal observables for and 
f 2fM are plotted as a function of the cut parameter y in Fig. 3. 





C22 


C33 


2.64 • 10~ 3 


5.61 • 10~ 4 


1.17- 10~ 8 





L 






fin 


1 


0.000 


0.000 




0.000 


1 


-0.023 


4 


0.000 


-0.023 


1 



Table 2 The variances ca = (OiOi) sm and the rescaled covariance matrix 
(OiOj)sM(cuCjj)^ 1 ^ 2 of the optimal observables for the measurement of the 

parameters /i ;2(U and d 7 " . 





fin 


flix 




optimal observable 


0.195 


0.423 


92.5 


T 


0.238 






V 




0.523 


101 



Table 3 1 s.d. accuracies Sg (32) with which the couplings g = / 2jtt , and 

can be measured assuming a sample of 10000 decays Z — > yU + /i~7 within 
the cut (22). Missing entries correspond to a vanishing sensitivity (5g = 00). 



We note that a search for CP-odd effects in the two-body decay Z — > /i+yU" 
is in practice not possible [2] since the direction of flight of the muons carries 
no CP-odd information, and the spin analyzing muon decays usually escape 
detection. Nevertheless, this reaction should be useful to study possible detec- 
tor effects which may feign non-zero mean values or asymmetric distributions 
of the correlations (17). 
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5 Z decays to r + r 7 

Compared to the muon case, CP tests in Z decays to t + t~7 follow a somewhat 
different line: The r momenta can not (yet) be measured directly. Instead we 
use the momenta of the charged decay products of the r leptons to build the 
CP-odd correlations T, V and their counterparts T*, V* which use the momenta 
in the t + t~ rest frame. As explained in Sect. 3 we assume, however, knowledge 
of the momenta for the optimal observables. 

In order to cover a large fraction (about 35%) of the r + r~ events, we in- 
vestigate all combinations of the following one-prong decay channels: 

(33) 
(34) 

(£ = e,/i). (35) 

For these decays the SM matrix elements are assumed. 

In Tables 4,5 we collect the upper limits on the CP-odd couplings which can 
be obtained by measuring the optimal observables and the correlations T, V, 
T*, and V*. The numbers N of events for the various channels correspond 
roughly to a total number of 10 4 Z — > t + t~7 events within the cut (22). 
Here we always assume that only one parameter g = dj,d^ , fi T , fn is different 
from zero and quote Sg calculated according to (31, 32) for the observable 
indicated. The measurement of the dipole moments relies essentially on the 
analysis of the r spins. The momentum distribution of the final state particles 
depends on the r polarization. For instance, the n + in the reaction (33) is 
emitted preferentially opposite to the r + spin direction, if we consider the 
decay in the r + rest system. The momentum correlations T,V,T*,V* which 
we use project on the CP-asymmetric part of the spin-spin-correlation of the 
r pair. Therefore they are rather insensitive to the electric dipole coupling dj 
which affects the polarization of only one r lepton (see Fig. lc). The weak 
dipole coupling d%, however, can be measured with an accuracy of the order of 
10~ 17 ecm assuming that 10 4 Z — > t + t~7 decays are available. The obvious gap 
between the sensitivities of the optimal observables and the simple correlations 
originates from the fact that the optimal observables as we calculate them 
work with the knowledge of the r momentum. At least for the measurement 
of the electric dipole moment this seems to be a basic ingredient. Since a 
measurement of the EDM and WDM is based on the r-spin information, we 
observe a rather strong dependence of the sensitivities in Tab. 4 on the r-decay 
channels, corresponding to their spin-analyzing qualities. 



-> vvt{k ± ) 
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The parameters f\ T , fir can be determined with an accuracy of 1 s.d. when 
they are of the order of about 0.5 (Tab. 5), given the assumed number of events. 
Here it should be advantageous to use the observables T and V which get only 
tiny contributions from the dipole moments of the tau lepton. These results are 
comparable to the muon case (Tab. 3). It is worth noting that the sensitivities 
are more or less independent of the tau decay channel: the measurement of 
fir, far is not substantially based on the spin analyzing qualities of the different 
decay modes. 









channel 


N 


opt. 


T 


V 




v* 


opt. 


T 


V 




V* 


7T — 7T 


100 


0.40 


91 


59 


8.4 


5.0 


0.054 


0.83 


8.3 


0.077 


0.67 


71 — p 


400 


0.41 


67 


41 


13 


3.2 


0.038 


0.63 


38 


0.051 


1.8 


p-p 


400 


0.85 


150 


83 


42 


5.9 


0.070 


1.3 


28 


0.092 


44 


l-l 


800 


0.41 


74 


51 


14 


4.8 


0.044 


0.73 


12 


0.063 


0.63 


e-7i 


600 


0.26 


120 


95 


7.5 


17 


0.030 


0.98 


3.4 


0.140 


0.46 


i-p 


1200 


1.5 


1400 


810 


53 


40 


0.13 


38 


23 


0.91 


4.3 


combined 


0.17 


40 


26 


4.8 


2.2 


0.018 


0.37 


3.0 


0.032 


0.32 



Table 4 1 s.d. accuracies with which <P T and d% can be measured 
for a given number N of events in various r decay modes 
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fir 


At 


channel 


N 


opt. 


T 


V 




v* 


opt. 


T 


V 


T* 


v* 


7T — 7T 


100 


0.61 


1.7 


7.1 


104 


45 


1.1 


3.4 


3.7 


36 


31 


7T — p 


400 


0.53 


0.85 


5.9 


59 


32 


0.88 


2.9 


1.9 


38 


36 


p- P 


400 


0.70 


0.88 


17 


60 


57 


1.3 


8.6 


2.0 


46 


68 


i-l 


800 


0.48 


0.60 


4.8 


29 


62 


0.87 


2.2 


1.4 


22 


34 


£-71 


600 


0.39 


0.66 


7.2 


20 


700 


0.68 


3.2 


1.5 


23 


140 


t-p 


1200 


1.4 


1.7 


130 


74 


340 


2.5 


64 


3.8 


190 


1200 


combined 


0.23 


0.34 


3.0 


15 


22 


0.4 


1.4 


0.78 


13 


18 



Table 5 1 s.d. accuracies with which f lT and fi T can be measured 
for a given number N of events in various r decay modes 



6 Conclusions 

In this article we have investigated the radiative decays Z — > £ + £~'~f (£ = p, r) 
with respect to the possibility of performing tests of the CP symmetry. The 
relevant CP-odd parameters are the EDM and WDM coupling parameters of 
the leptons, dj, df and 4-point coupling parameters fa. For the muon case 
the EDM does not contribute to CP-odd correlations in Z decays to leading 
order. But in the effective Lagrangian approach eft can be bounded by the 
limits obtained in the direct search for an EDM of the muon. On the other 
hand, the radiative Z decays offer the possibility to measure the muon's WDM 
and its 4-point couplings / 2(U . The accuracies obtainable with a realistic 
number of events are collected in Tab. 3. 

For the r-lepton one can in principle measure all four CP-odd parameters 
dl, , fir, far using radiative Z decays. The measurement of d% using about 
10 4 Z — > t + t~7 decays has a 1 s.d. sensitivity (cf. Tab. 4) 

5d z T = 0.018. (36) 

This is, as expected, weaker than the bound (15) obtained from Z — > t + t~ 
decays. Radiative Z decays offer the possibility to perform a true measurement 
of the EDM coupling parameter dj. To date no direct determination of dl 
exists. Another reaction where dj can be measured is e + e~ — > t + t~ away from 
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the Z resonance (cf. [7]). With 10 6 r-pairs at y/s = 10 GeV for instance one 
gets a sensitivity: 



5<F T = 0.043. (37) 

This is to be compared to 5dJ = 0.17 or 2.2 in Tab. 4 using the optimal 
observable or V*, respectively. Finally, the parameters fi T , for are special to 
Z — > t + t~7 decays and can be measured with accuracies of order 0.5 with the 
number of events given in Tab. 5. 

An investigation of radiative Z decays can thus give information of CP-odd 
coupling parameters of leptons which is complementary to the one obtainable 
from other sources. We have found that for the observables and cuts we con- 
sidered correlations between the CP-odd coupling parameters are negligible. 
Thus we investigated each coupling parameter separately. In an analysis of 
real experimental data other cuts may be used and smearing due to detec- 
tor effects will play a role. One will again be confronted with the question 
of correlations and in a sophisticated analysis one will have to consider the 
parameters fu, fee, cQ, df simultaneously for each lepton t. It should then be 
advantageous to use couplings and observables which give diagonal covariance 
matrices. The general methods to achieve this which are explained in [27] can 
easily be adapted for our case here. 

Independently of CP-violation it should be interesting to study the reaction 
Z — > £ + £~'~f in the phase space region where the photon is "hard". Any new 
effects in lepton physics, like substructures or excited leptons, could lead to 
contact interactions between Z,£ and 7. These should best show up in "hard" 
3-body decays, i.e. in regions of phase space where the photon has high energy 
and is well separated from £ + and £~. There the SM contribution is relatively 
small. An investigation of the ratio 

= r s „ ( z FF% < 38 > 

as function of the cut parameter y (22) should already be revealing. As an 
example we have calculated this ratio for our ansatz SM plus Cqp (2) assuming 
that only /v£ an d f\e are different from zero. To a very good approximation 
(cf.[9]) these couplings add to T(Z — > £ + £~) a term ATcp(fve, fAf,y) propor- 
tional to /v£ + f\t- Thus we plot 

R{y) - 1 = Ar C p(/ w ,/ A £;2/) , 3g ^ 

Me + fh IWZ-* €+€-7) + 

in Fig. 4 (solid line). We see that a value fy e + f\ t = 1 will lead to deviations 
of R(y) from 1 of up to about 10% at y = 0.2. An investigation of 
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plotted in Fig. 4 (dashed line) shows that the absolute contribution of CP- 
odd parameters to the width decreases rapidly for increasing values of the cut 
parameter y. Here ATcp(fve, fkt) is the addition to the width for cut parameter 
y = as given in (10). 

The calculations reported in this article required a considerable amount 
of symbolic and numerical computation. We used the new programming lan- 
guage M [28] for the symbolic manipulations including the expansion of traces of 
Dirac-matrices, simplification of the resulting large expressions, and generation 
of optimized Fortran program code. The numerical phase space integrations 
(with dimensions up to 11) were done with the Fortran routine VEGAS. The 
programs which we used as well as Monte Carlo event generators for the reac- 
tions discussed in this article can be obtained from the authors (World Wide 
Web address: http://www.thphys.uni-heidelberg/~ overmann). 
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Figure Captions 



Figure 1 

Feynman diagrams contributing to Z — > 

a) One of the lowest order Standard Model diagrams 

b) Contribution of f^e and fye 

c) Contribution of the electric dipole moment dj 

d) Contribution of the weak dipole moment df 

Figure 2 

Contour plot of the 2 s.d. errors on / V/1 , /a^ obtainable ideally from the mea- 
surement of the correlations T, V from Tab. 3 (stripes) and from the contribu- 
tion to the width (circle) (c.f. (11,14)) assuming 10000 fi + [i~'y events and zero 
mean values. 

Figure 3 

1 s.d. accuracies for the determination of fi^ and / 2/J for different values of 
the phase space cut parameter y. Here we assume a number of /i + /z~7 events 
corresponding to 10 4 for y = 0.03. 

Figure 4 

The ratio (R(y) — l)/(/v* + fie) ( s °lid li ne ) as given in (38,39), and the ratio 
r(y) (dashed line) defined in (40) as function of the cut parameter y (22). 
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